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Abstract. Guided by non-relativistic Effective Field Theory (EFT) we classify the most 
general spin-dependent interactions between a fermionic Weakly Interacting Massive Par¬ 
ticle (WIMP) and nuclei, and within this class of models we discuss the viability of an 
interpretation of the DAMA modulation result in terms of a signal from WIMP elastic scat¬ 
terings using a halo-independent approach. We find that, although several relativistic EFT’s 
can lead to a spin-dependent cross section, in some cases with an explicit, non-negligible 
dependence on the WIMP incoming velocity, three main scenarios can be singled out in the 
non-relativistic limit which approximately encompass them all, and that only differ by their 
dependence on the transferred momentum. For two of them compatibility between DAMA 
and other constraints is possible for a WIMP mass below 30 GeV, but only for a WIMP veloc¬ 
ity distribution in the halo of our Galaxy which departs from a Maxwellian. This is achieved 
by combining a suppression of the WIMP effective coupling to neutrons (to evade constraints 
from xenon and germanium detectors) to an explicit quadratic or quartic dependence of the 
cross section on the transferred momentum (that leads to a relative enhancement of the 
expected rate off sodium in DAMA compared to that off fluorine in droplet detectors and 
bubble chambers). For larger WIMP masses the same scenarios are excluded by scatterings 
off iodine in GOUPP. 
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1 Introduction 

Many underground experiments are currently searching for Weakly Interacting Massive Par¬ 
ticles (WIMPs), which are the most popular candidates to provide the Dark Matter (DM) 
which is believed to make up 27% of the total mass density of the Universe[l], and more 
than 90% of the halo of our Galaxy. One of them (DAMA[2]) has been observing for more 
than 15 years a yearly modulation effect in the low part of its energy spectrum which is 
consistent with that expected due to the Earth rotation around the Sun from the elastic 
scattering of WIMPs off the sodium iodide nuclei that constitute the crystals of its scin¬ 
tillators. Many experimental collaborations using nuclear targets different from Nal and 
various background-subtraction techniques to look for WIMP-elastic scattering (LUX [3], 
XENON100[4], XENON10[5], KIMS[6, 7], CDMS-Ge[8], CDMSlite [9], SuperCDMS[10], 
SIMPLE[11], COUPP[12], PICASSO[13], PICO-2L[14]) have failed to observe any anomaly 
so far, implying severe constraints on the most popular WIMP scenarios used to explain the 
DAMA excess. This is particularly compelling when the KIMS bound [6] (which is obtained 
using Csl crystals) is taken at face value, and when the WIMP mass is large enough to rule 
out the possibility that the DAMA effect is mainly due to WIMP scatterings on sodium. 
In that case both DAMA and KIMS should observe the same WIMP-iodine process and no 
room is left to model building to reconcile the discrepancy (although several experimental 
uncertainties might still be advocated to question the robustness of the KIMS bound^). How¬ 
ever, when the WIMP mass is small enough it is possible to assume that the contribution of 

^The raw background of KIMS is about a factor of three larger than that of DAMA, so that KIMS sensitivity 
relies on background subtraction techniques that can potentially introduce systematic errors. Moreover KIMS 
has recently published a new measurement of its quenching factor [15] significantly smaller than that used for 
its published bounds[6, 7] and reducing its sensitivity to WIMPs scatterings at low energy. 
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WIMP-iodine scatterings in the expected event rates of both DAMA and KIMS vanishes and 
that the DAMA effect is due to scatterings on sodium. In this case the constraints on DAMA 
can be potentially relaxed by considering models where the expected WIMP scattering rate 
on sodium is enhanced compared to that off the nuclei used to obtain the experimental 
bounds. 

Moreover, the energy dependence of the WIMP-induced scattering spectrum depends 
on the velocity distribution f(v) of the incoming WIMPs. Historically, for the latter the 
isothermal sphere model has been assumed, i.e. a thermalized non-relativistic gas described 
by a Maxwellian distribution whose temperature Vrms — 270 km/sec is determined from the 
galactic rotational velocity by assuming equilibrium between gravitational attraction and 
WIMP pressure. However, while the presence of a thermalized component in the WIMP 
halo of or Galaxy is conhrmed by numerical simulations the detailed merger history of the 
Milky Way is not known, allowing for the possibility of the presence of sizeable non-thermal 
components for which the density, direction and speed of WIMPs are hard to predict. 

In light of the situation summarized above several new directions have been explored 
in the recent past both to remove as much as possible the dependence on specific theoretical 
assumptions (both of particle-physics and astrophysical origin) from the analysis of DM 
direct detection data and to extend its scope to a wider class of models. Starting from 
[16] a general strategy has been developed [17-19] to factor out the dependence on f{v) of 
the expected WIMP-nucleus differential rate dR/dE^ at the given recoil energy Er. This 
approach exploits the fact that dR/dEn depends on f{v) only through the minimal speed 
Vmin that the WIMP must have to deposit at least Er, i.e.: 


dR 

dEji 




/ M 

J \v\^V^in 1^1 


( 1 . 1 ) 


By mapping recoil energies Er into same ranges of Vmin the dependence on ij{vmin) and so on 
f{v) cancels out in the ratio of expected rates on different targets. Since the mapping between 
Er and Vmin depends on the nuclear mass the factorization of rj{vmin) is only possible in the 
case of detectors with a single nuclear target, or for which the expected rate is dominated by 
scatterings on a single target. However in the following we will extend this procedure to check 
the compatibility between a candidate signal (such as that from the DAMA experiment) and 
a null experiment also in the case when the latter contains different targets. 

On the particle-physics side one of the most popular scenarios for WIMP-nucleus scat¬ 
tering is a spin-dependent interaction where the WIMP particle x is a fermion (either Dirac 
or Majorana) that recoils on the target nucleus T through it coupling to the spin S'tv of 
nucleons N = {p,n): 


Cint <xS^-SN = C^S^-Sp + C^S^-Sn. ( 1 - 2 ) 

The above effective Lagrangian can arise as the low-energy limit of a relativistic axial 
coupling (the incoming WIMP velocity in the halo of our Galaxy is ex¬ 

pected to be of order 10“^c). One of the main motivations of such scenario is the fact that 
the most stringent bounds on the interpretation of the DAMA effect in terms of WIMP- 
nuclei scatterings arise today from detectors using xenon (LUX[3], XENON100[4]) and ger¬ 
manium (SuperCDMS[I0]) whose spin is mostly originated by an unpaired neutron, while 
both sodium and iodine in DAMA have an unpaired proton: as already pointed out by 
some authors [20, 21], if the WIMP effective coupling to neutrons c"" is suppressed com¬ 
pared to that on protons d’ this class of bounds can be evaded. However this scenario is 
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presently constrained by droplet detectors (SIMPLE[11], COUPP[12]) and bubble chambers 
(PICASSO [13], PICO-2L[14]) which all use nuclear targets with an unpaired proton (in par¬ 
ticular, they all contain while SIMPLE contains also and and COUPP uses 
also As a consequence, this class of experiments have been shown to rule out the 

scenario of Eq. (1.2) also for <C when standard assumptions are made on the WIMP 
local density and velocity distribution in our Galaxy [14, 21]. 

The spin-dependent scenario has been further extended in Refs. [22, 23], where a 
generalization of the coupling of Eq.(1.2) has been discussed. In this scenario the WIMP 
particle couples to nucleons through a relativistic axial coupling through the 

exchange of a light pseudoscalar, leading to a non-relativistic coupling {Sy. ■ q){S]y ■ q) which 
depends on the projection of the nuclear spin component along the transferred momentum q. 
This scenario, which combines a spin-dependent nuclear response function mainly sensitive 
to unpaired nucleons (leading to relaxed constraints from Ge and Xe detectors on DAM A 
when c” <C c^) to a non-standard dependence on the transferred momentum q, has also been 
shown to be in tension with bubble chambers and droplet detectors, albeit to a lesser extent, 
at least for standard assumptions on the WIMP velocity distribution [21]. 

In the present paper we wish to elaborate further on the possibility that the non- 
relativistic effective coupling between WIMPS and nuclei is due to a combination of a spin- 
dependent nuclear response function and a generalized momentum dependence. In particular, 
we wish to assess the viability of this kind of scenario by departing from the Isothermal Sphere 
model and adopting the previously mentioned halo-independent approach. 

Our main motivation is the simple observation that models such as that of Refs. [22, 23], 
where the WIMP-nucleus cross section depends on additional factors g” (with n >0) of 
the transferred momentum, tend to suppress the sensitivity to the DAMA effect of droplet 
detectors and bubble chambers because, due to a combination of low energy thresholds and 
light target masses, a critical fraction of the WIMP-scattering events expected in the latter 
have transferred momenta below those related to the DAMA effect. For this reason, scenarios 
that combine a spin-dependent coupling mostly to protons to additional factors g” in the cross 
section are expected to have better chances to reconcile the DAMA effect to the constraints 
from all other detectors. This will be confirmed in our analysis. 

In order to develop this program in a systematic way, we will rely on the analysis 
of Ref. [24], where the most general non-relativistic Hamiltonian H describing the elastic 
scattering of a fermionic WIMP particle off nuclei was written in terms of the sum of all 
the terms invariant by Galilean transformations. Using this approach, we will single out all 
the possible spin-dependent interactions compatible to Galilean invariance that will be the 
subject of our phenomenological analysis. 

Our paper is organized as follows. In Section 2 we use the non-relativistic Effective 
Field Theory approach of Ref. [24] to single out the generalized spin-dependent interactions 
we wish to discuss in the rest of our paper, and set our notations; in Section 3 we provide 
the formulas to calculate expected rates for WIMP-nucleus scattering in our scenarios; in 
Section 4 we summarize the main ingredients for the halo-independent analysis; in Section 
5 we introduce a compatibility factor between DAMA and other experiments and provide a 
recipe to extend it to the case of constraints involving more than one nuclear target. The 
results of our quantitative analysis are given in Section 6, while Section 7 is devoted to our 
conclusions. Finally, in the Appendices we provide some technical details of our procedure 
and summarize the experimental inputs used in the analysis. 
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2 Generalized spin—dependent interactions 

The most general WIMP-nucleus spin-dependent interactions can be single out by making 
use of the non-relativistic EFT approach of Ref. [24]. According to that analysis the most 
general Hamiltonian density for the process can be expressed in terms of a combination of 
the following five Hermitian operators, which act on the two-particle Hilbert space spanned 
by tensor products of WIMP and nucleon states: 

l^AT iq v-^ Sn , (2.1) 

where is the identity operator, g is the transferred momentum, and Sn are the WIMP 
and nucleon spins, respectively, while v-^ = v + ^ (with the WIMP-nucleon reduced 

mass) is the relative transverse velocity operator satisfying g = 0. Including terms that 

are at most linear m. Sn, Sy_ and v'^ only the following 15 non-relativistic quantum mechanical 
operators can be constructed out of (2.1): 
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ruN mN 

so that the the most general Hamiltonian density describing the WIMP-nucleus interaction 
can be written as: 

15 

= E (2.3) 

r=0,l k=l 


where = 1, denote the the 2x2 identity and third Pauli matrix in isospin space, 

respectively, and the isoscalar and isovector (dimension -2) coupling constants and c^, are 
related to those to protons and neutrons and by = (c° -|- c^)/2 and = {c^ — c\)/2. 

The basic assumption of Ref. [24] is that the nuclear interaction is the sum of the inter¬ 
actions of the WIMPs with the individual nucleons in the nucleus, and so neglects any effect 
involving two or more nucleons. Notice that while such effects are usually expected to be 
small, they can become non-negligible in particular circumstances, such as when the leading 
terms is suppressed by a cancellation between the and couplings [25, 26]. With this 
simplification the WIMP scattering amplitude on the target nucleus T can be written in the 
compact form: 
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In the above expression and Jt are the WIMP and the target nucleus spins, respectively, 
q = \^ while the ’s are WIMP response functions (that we report for completeness in 
Eq.(A.l)) which depend on the couplings as well as the transferred momentum q and on 
the square of the transverse velocity in the reference frame of the center of mass of the target 
nucleus which can be written as: 


= 4 - u 


2 

min' 


(2.5) 


In particular, the quantity: 


■ = 


r 


rriTEji 


( 2 . 6 ) 


2^t 

represents the minimal incoming WIMP speed required to impart the nuclear recoil energy 
Eji, while vt = lilrl is the WIMP speed in the reference frame of the nuclear center of 
mass, rriT the nuclear mass and jix the WIMP-nucleus reduced mass. In equation (2.4) 
the (9^)’s are nuclear response functions and the index k represents different effective 
nuclear operators, which, crucially, under the assumption that the nuclear ground state is 
an approximate eingenstate of P and CP, can be at most eight: following the notation 
in [24], k=M, <h"M, S", S', A,AS'. For the target nuclei T used in most direct 

detection experiments the functions (g^), calculated using nuclear shell models, have 
been provided in Refs. [24, 27]. 

In particular, the usual spin-dependent nuclear scattering process corresponds to R = 
X)r=o 1 The corresponding form factor is usually written in the literature in terms 

of the nuclear spin structure functions as [28]: 


S{q^) = (clfSooiq^) + 4 cISoi{q^) + (clfSooiq^). (2.7) 

At zero momentum transfer this form factor is normalized to a combination of Jt and of the 
expectation values < Sp^n >=< T|S'p^„|T > of the projection along the 2 axis of the total 
spin operators Sn,p = X) ^n,p of protons and neutrons on the nuclear state |T >: 

^ ^ 1 (2jT + l)(jr + l) ^^ ^ ^ ^ (2.8) 

TT jT 

In terms of the functions of Ref. [24] the structure functions Stt-' are given by: 


Soo{q^)=Wi?,iq^) + W^?{q^) 

Suiq^) =Wll{q^) + Wi,}{q^) 

Soiiq^) = 2 [Wil{q^) + Wl}{q^)] . (2.9) 

The above decomposition of the spin-dependent structure functions in terms of the two dif¬ 
ferent non-relativistic nuclear response functions S" and S' traces back to the decomposition: 

S^- Sn = {Sx ■ q/q){SN ■ q/q) + {S^ x q/q) ■ {Sn x q/q). ( 2 . 10 ) 

In particular, S" corresponds to the component of the nucleon spin along the direction of 
the transferred momentum q while S' to that perpendicular to it (for this reason one has 
W^T'{q^) ~ 2W^T!{q^) when q^ ^ 0). 
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Table 1. Nuclear response functions corresponding to each coupling, for the velocity-independent 
and the velocity-dependent components parts of the WIMP response function, decomposed as in 
Eq.(2.11). In parenthesis the power of q in the WIMP response function. 


Among the ’s, which represent the most general nuclear response functions for 
WIMP-nucleus scattering, only and W£T have the property of being suppressed for 

an even number of nucleons in the target, allowing to reconcile the DAMA effect to xenon 
and germanium detectors by assuming c” —)-0 (all other structure functions vanish in each 
different target for a different ratio r* = c”/c^). The correspondence between models and 
nuclear response functions can be directly read off from the WIMP response functions 
and (see Eq.A.l). In particular, using the decomposition: 
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( 2 . 11 ) 


this correspondence is summarized in Table 1. We see from this table that the models 
which couple to the nucleus either through W^T {q^) or W^, {q^) are Oi for i=3, 4, 6, 7, 
9, 10, 12, 13, 14, 15. However, for models Oi with i=3, 12, 13, 15 the spin-dependent 
coupling is velocity-suppressed, while the dominant velocity-independent term couples to 
another nuclear response function. This implies that the truly “spin-dependent” couplings 
are only Oi with i=4, 6, 7, 9, 10, 14. In Table 2, which is adopted from [24], we show which 
class of relativistic Effective Field theories can lead to them as a non-relativistic limit. For 
each relativistic model, in the last column we provide the dependence of the WIMP-nucleus 
scattering cross section on g, {q^) and W^J, {q^)- 

Irrespective of their specific velocity distribution, WIMP velocities are expected to be 
of order 10~^, so in line 4 of Table 2 the velocity-dependent term is naturally suppressed. 
For this reason we have neglected it. However, notice that in line 8 the coupling term 
leads to a velocity-suppressed cross section oc (v^)^ without the cor¬ 
responding velocity-independent contribution. The same thing can be achieved by appropri- 
ate linear combinations of the operator either with mwiMP ^ 

Xi^iiv mwiMP ^ as shown in lines 2 and 3 of the same Table. In Table 1 we ordered 

the models with powers of q^: actually, as we will see in the following, the scaling with q^ 
will be the most relevant property to interpret our numerical results. 

With the exception of line 4 and line 11, all the models of Table 2 correspond to one 
of the non-relativistic quantum-mechanical operators of Eq.(2.2). This also holds for the 
model of line 4 that can be well approximated by Og. So in the following we will refer to 
each of these scenarios with the corresponding Oi. On the other hand, the model of line 
11 is given by a superposition of O4 and 0% with both contributions of the same order. In 


- 6 - 














Relativistic EFT 

Nonrelativistic limit 

E,Oi 

cross section scaling 

1 


—• Sj\; 

- 4 O 4 

w-t ( 9 b + w-r' ( 9 b 

2 

‘^Xl'^xN-iij.p N+ 

—45"AT • 

- 4 O 7 


3 

— t;;— -xiV7^75iV’ 

—4iS/v • 

- 4 O 7 

(«4)2wjr'(9b 

4 


— 25^;- • i7^ + 

+ — - - - iS-x • (S]<^ X 

2 O 7 + 2 O 9 

m-wiMP 

~ 2 O 9 

'^WTMP 


5 

-i—X-^7^ 75 


4 

9^W-7'(9b 

6 

XX**75X Vio-JV 

■ C ^ ^n) 

_4j21iV Og 

mj^ 


7 


i 

m/V 

Oio 

q'^W^r/iq^) 

8 

iXio-ij.v 75 X V7,i75 iV 

■ ^x)('"t ■ Sw) 

-4^0i4 

(v^Pq^W^f (9b 

9 

X75XV7^JV 

^ ■ Sv U ■ SiV 

^WTMP 

^ Os 

^WTMP ° 

9^W-7/(9b 

10 


44ySx4ySU 

4^06 

'^M 

9^W-,V(9b 

11 


^ ®x) ■ (^ X Sn) 

4 ( -4-04 - ^Og] 

9^W-7'(9b 


Table 2. Relativistic Effective Field Theories for a Dark Matter fermionic WIMP x having as a 
low-energy limit a generalized spin-dependent x-nucleus elastic scattering. In all these scenarios the 
scattering cross section vanishes for all nuclei with an odd number of protons if =0 and for all nuclei 
with an odd number of neutrons if cf=0. Some of the interaction terms in the second column contain 
an arbitrary scale tom to ensure correct dimensionality. 


Section 6 we will conventionally refer to this model as 046- In summary, we conclude that the 
most general spin-dependent models relevant from the point of view of the phenomenology 
of direct detection are seven: O 4 , O7, O9, Oio, O 14 , Oq, O 4 Q. We will quantitatively discuss 
their phenomenology in Section 6. 


3 WIMP direct detection rate 


For a given recoil energy imparted to the target the differential rate for the WIMP-nucleus 
scattering process is given by: 
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(3.1) 


where pwiMP is the local WIMP mass density in the neighborhood of the Sun, Nt the 
number of the nuclear targets of species T in the detector (the sum over T applies in the case 
of more than one nuclear isotope), /(hr) the WIMP velocity distribution, while 
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with the squared amplitude in parenthesis given explicitly in Eq.(2.4). Using Eqs (2.11,3.1,3.2) 
one gets: 
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In Eq.(3.3) we have factorized a reference cross section aref which in the case of no momentum 
dependence of the scattering amplitude can be identified with the large-distance total cross 
section for a point-like nucleus. However in the case of scattering amplitudes which vanish 
with the quantity Uref is just a conventional factor with dimension -2 which must be 
the same for all nuclei in order to cancel out in the comparison of the expected rates in 
different detectors. For all the interactions listed in Table 2 only one of the couplings is 
different from zero. In the following we choose to factorize for each case the corresponding 
conventional cross section cr^e/ = 

Eq. (3.3) implies that, among the models summarized in Table 2, those in columns 2, 
3 and 8 have a direct detection rate proportional to ^{vmin)-, while all the others scale as 
fj(vmin)- In both cases the same function can be factorized in all nuclei, allowing to adopt 
the halo-independent procedure outlined in the next Section. 


4 Halo—independent analysis 

We summarize in this Section the main formulas used in our subsequent analysis to factorize 
the halo-dependence from direct-detection data. 

In a realistic experiment the recoil energy Er is obtained by measuring a related detected 
energy E' obtained by calibrating the detector with mono-energetic photons with known en¬ 
ergy. However the detector response to photons can be significantly different compared to the 
same quantity for nuclear recoils. For a given calibrating photon energy the mean measured 
value of E' is usually referred to as the electron-equivalent energy E^e and measured in keVee. 
On the other hand Eji (that represents the signal that would be measured if the same amount 
of energy were deposited by a nuclear recoil instead of a photon) is measured in keVnr. The 
two quantities are related by a quenching factor Q according to Eee = Q{Er)Et- More¬ 
over the measured E' is smeared out compared to Egg by the energy resolution (a Gaussian 
smearing Gauss{Ef,e\E', arms{E')) = 1/{armsV^)exp[—{E^ — Ege)"^/{‘icr^ms)] with standard 
deviation arms{E') related to the Full Width Half Maximum (FWHM) of the calibration 
peaks at E' by FHWM = 2.35arms is usually assumed) and experimental count rates de¬ 
pend also on the counting efficiency or cut acceptance e(E'). Overall, the expected differential 
event rate is given by: 


JD roo -j TD 

— = 6(E') dEe,GaussiE,e\E',arms{E'))-^^^^, (4.1) 

with dR/dEji given by Eq.(3.1). Then the number of events in a given interval of the detected 
energy E' is simply: 
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Combining Eqs.(4.1,4.2) with Eq. (3.3) and changing variable from Er to Vmin the expected 
rate can be cast in the compact form: 
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(4.5) 


In the present analysis we will assume that the yearly modulation effect observed by 
the DAMA experiment is due to the time-dependence of r) or ^ due to the rotation of the 
Earth around the Sun. So we also introduce the modulated amplitudes: 


Pliymin) — p(,PminR — tmax') ViPniinR — tmin)i 
Cl(^mm) = ^{,'^minR — tmax) ^{p^minit — 


(4.6) 


with tyyiax a-i^d tmm the times of the year corresponding to the maximum and to the minimum 
of the signal. The modulated amplitudes measured by DAMA that we will use in Section 
5 to get estimations of fji and were obtained by fitting the data with a time dependence 
parametrized with the functional form: 


R 


[E(,Eb — -^[EbEb + -^^EbEb COs[a;(t - tmax)], 


(4.7) 


which is expected when the isothermal sphere model is adopted for the WIMP speed dis¬ 
tribution. Moreover, according to the DAMA analysis, the modulation phase is in (rough) 
agreement to what is expected for the isothermal sphere model, i.e. tmax — 2 June (and, 
consequently, tmm — 2 December). However, a truly halo independent analysis would require 
to reanalyze the experimental data allowing for functional forms of general nature beyond 
Eq.(4.7)^. Nevertheless, in the following to estimate the fji, functions we will make directly 
use of the published DAMA modulation amplitudes, implicitly assuming that indeed those 
estimations do not differ significantly from what one would obtain using a more general form 
to analyze the data. With this important proviso, in the present analysis we wish to make 
the smallest possible number of assumptions on the functions r), They turn out to be the 
same, and in the case of i) they reduce to: 


^For speed distributions beyond the Isothermal Sphere the time dependence may differ from that of 
Eq.(4.7). For instance, it has been shown that in anisotropic extensions of the isothermal sphere the phase 
and even the period may vary [29]. 
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The first condition descends from the definitions (3.4), that implies that both fi[vmin) 
and i{vmin) are decreasing functions of Vmin- The second is an obvious consequence of the 
fact that r)i, are the modulated parts of fj and The last condition reflects the requirement 
that the WIMPs are gravitationally bound to our Galaxy. In the following we will assume 
that the WIMP halo is at rest in the Galactic rest frame and we will adopt as the maximal 
velocity of WIMPs Uesc=782 km/s in the reference frame of the laboratory, by combining the 
reference value of the escape velocity km/s in the galactic rest frame with the 

velocity ^0=232 km/s of the Solar system with respect to the WIMP halo. 

So, given an experiment with detected count rate A^exp in the energy interval E'^ <E' < 
E 2 the combination: 


< fj >= 


f-dv. 


I 


lo dv min '^[E[,E^]iVn 


N, 


exp 


lo dv min T^[E[,E'2]{Vn 


(4.9) 


can be interpreted as an average of the function fj{vmin) in an interval Vmin,i < Vmin < Vmin, 2 , 
with an analogous definition for the average of the function The Vmin interval is 

defined as the one where the response function TZ is “sizeably” different from zero (we will 
conventionally take the interval Vmin[ER{Eee,i)] < Vmin < Vmin[ER{Eee, 2 )] with Eee,i = 
E'l — (JrmsiE'i), Eee ,2 = (-^2)> i-®- the E' interval enlarged by the energy resolution). 


5 Constraints from direct detection experiments and compatibility factor 

In this Section we will give explicit formulas for the case of a velocity-independent cross 
section where the expected rate is given by Eq.(4.3) and depends on fj{vmin)- Clearly, the 
velocity-dependent case of Eq.(4.4) with dependence on ^{vmin) is analogous. 

Eollowing the halo-independent procedure outlined in Section 4 it is straightforward, for 
a given choice of the DM parameters, to obtain estimations of the modulation- 

amplitude function fji{vmin) averaged in appropriately chosen Vmin,i intervals mapped from 
the DAMA experimental annual modulation amplitudes. Using the condition fji{vmin) < 
v{vmin) this allows to get lower bounds on the fj(vmin) function, which can be compared to 
upper bounds fij^nm on the same quantity derived from the data of the experiments that have 
reported null results. Among the latter, the KIMS experiment uses iodine, one of the same 
targets of DAMA {Nal for DAMA vs. Csl for KIMS). KIMS has reported an upper bound 
on the unmodulated WIMP rate in the low-part of its measured spectrum that is about a 
factor of two smaller than the modulated amplitude measured by DAMA in the same energy 
interval [6] . This implies that, as far as the WIMP rate in both experiments is dominated 
by scatterings on iodine, the two results cannot be reconciled by changing either the halo 
function or the scaling law of the cross section. In this case the only way to reconcile the two 
conflicting results is to look deeper in the possible sources of systematic errors, including the 
many uncertainties connected to quenching factors, atomic form factors, background cuts 
efficiencies, etc. Moreover, as already pointed out, as far as DAMA is concerned the halo- 
function factorization is only possible if WIMPs scatter predominantly on a single target 
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nucleus. For both reasons in the following we will restrict our analysis to a range of the 
WIMP mass for which the DAMA effect can be consistently explained with scatterings on 
sodium targets only. In the popular case of an Isothermal Sphere for the velocity distribution 
this is achieved for mwiMP < 20 GeV. However this range for mwiMP can be larger if the 
WIMP velocity distribution is different from a Maxwellian and if the only assumptions on 
the halo function are those of Eq.(4.8). In particular it can be maximally extended by 
assuming that the function f] is the minimal one compatible to the DAMA signal. Assuming 
that the latter corresponds for sodium to the Vmin range < yDAMA,Na 

and for iodine to the range the most conservative assumption 

on fj is that it vanishes for Vmin > particular, this is consistent with the 

minimal requirement that the halo function r] is a decreasing function Vmin-)- Then, as 
long as > V 2 , one can conservatively assume that scatterings on iodine are 

below the energy threshold of DAMA and do not contribute to the modulation effect. Using 
Eq.(2.6) this corresponds to mwiMP < 60 GeV, assuming the energy range 2 keVee < E' < 
4 keVee for the DAMA modulation region, enlarged by an energy resolution ~ 0.448 x y/W 
keVee [2]. In the same range of mwiMP also scatterings on iodine and caesium are below 
the threshold of 3 keVee in KIMS (assuming the latest measurement of the KIMS quenching 
factor [15], which is similar to that of DAMA). As a consequence of this in the following we 
will restrict our analysis to mwiM P 60 GeV, assume only scatterings on sodium in DAMA 
and neglect the KIMS constraint. 

Quantitatively, for a given choice of the WIMP mass mwiMP, of the interaction term 
from Table 2 and of the corresponding ratio r*, the compatibility between DAMA and all the 
other results can be assessed introducing the following compatibility ratio[26]: 

/ ~DAMA-Na , 

E>{mwiMP,ri) = max ^- 5 - - , (5.1) 

iGsignal y mmj<j rjj^um J 

where ai represents the standard deviation on a.s estimated from the data, i G 

signal means that the maximum of the ratio in parenthesis is for Vmin,i corresponding to the 
DAMA excess, while, due to the fact that the function fj is non-decreasing in all velocity bins 
Vmin,i, the denominator contains the most constraining bound on fj for Vmin,j < Vmin,i- The 
latter minimum includes all available bounds from scintillators, ionizators and calorimeters 
(see Appendix B for a summary of the experimental inputs used in our analysis). Specihcally, 
compatibility between DAMA and the constraints included in the calculation of Eq.(5.1) is 
ensured if P < 1. Notice that, by combining different Vmin,i bins, the above definition takes 
into account the most general momentum dependencies predicted by the various effective 
interaction terms of Table 2. 

The procedure outlined above cannot include bubble chambers and droplet detec¬ 
tors which are only sensitive to the energy threshold, such as SIMPLE[11](using C 2 CIF 5 ), 
GOUPP[12] (using CF 3 I), PIGASSO[13](using C4F10) or PIGO-2L[14] (using CgFg), since 
in this case it is not possible to map the corresponding bounds to arbitrary velocity bins. 
Notice that these experiments contain all proton-odd elements (fluorine, chlorine and io¬ 
dine), so they are particularly relevant for our analysis since they are complementary to 
germanium and xenon detectors which are neutron-odd. Moreover, another complication in 
using the compatibility factor defined in Eq.(5.1) to assess agreement between the DAMA 
excess and this class of experiments is that they all contain different nuclear targets, so that 
it is in general not possible to factorize the fj function and it is not trivial (as for instance 
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for the case of a coherent interaction, where the cross section scales with the square of the 
atomic mass number) to single out particular kinematic situations where one of the targets 
dominates, this also because of the momentum-dependence of some of the interactions. So, 
in order to impose the bounds from this class of experiments in a halo-independent way, 
we need to generalize the compatibility factor dehnition of Eq.(5.1). In order to do so, we 
proceed in the following way: i) we use the experimental DAMA modulation-amplitudes to 
get a conservative piecewise estimation of the minimal fji modulated halo function 

compatible to the signal; ii) we obtain the corresponding estimation of the unmodulated part 
by requiring that it is a decreasing function of Vmin with fj^^\vmin) > 
iii) in compliance with (4.8) and with the goal of obtaining a conservative bound, we re¬ 
quire that the function fj is the minimal one able to explain the DAMA effect, so we assume 
> V 2^^^’^°‘)=0 (this is also consistent to dominance of scatterings on sodium in 
DAMA in the mass range mwiMP < 60 GeV, as explained above); iv) we then use 
to directly calculate for each experiment among /c=SIMPLE, COUPP, PICASSO and PICO- 
2L and for each energy threshold Eth,i the expected number of WIMP events a,nd 

compare it to the corresponding upper bound (gee Appendix B for further details). 

Then, a straightforward generalization of the compatibility factor of Eq.(5.1) is: 

( j^expected \ 

E>{mwiMP,ri), 1 ■ (5.2) 

Notice that the above procedure is very general, since it allows to get a halo-independent 
bound in the case of any experiment that does not observe a signal. In particular this means 
that the requirement that one target dominates the expected rate is only needed to get 
estimations of the halo function from those experiments that observe an excess, but is not 
indispensable to get a conservative bound for those reporting a null result. In the following 
Section we will use the above dehnition of the compatibility factor to explore the WIMP 
parameter space for the generalized spin-dependent interactions introduced in Section 2. 

6 Results 

Each of the seven phenomenological models Ot (with i=4, 6, 7, 9, 10, 14, 46) listed at the end 
of Section 2 depends on three parameters: the WIMP mass mwiMP and the two components 
(c^,c”) of the corresponding effective coupling. However, following the halo-independent 
approach summarized in Section 4, we get our estimations of the fj or ^ functions factorizing 
out the reference cross section aref = a consequence of this the coupling cf 

cancels out in our dehnition (5. 1,5. 2) of the compatibility factor D. So in each model the 
phenomenology depends only on two parameters: mwiMP and the ratio r* (in the case of 
model C >46 we factorize a^ef = Cg/ij^-Ar/Tr). 

Moreover, the pairs of models Oi-O-j and Og-Ou only differ for an additional (u^)^ 
factor in the cross section (see Table 2), which is absorbed in the dehnition of the halo 
function So in each pair the two models have identical compatibility factors D (in each 
case the estimations for the velocity-dependent model are simply rescaled by a common 
factor with respect to the corresponding fji in the velocity-independent one) . The bottom 
line is that only hve relevant cases remain: O 4 , Oq, Og, Oio and O 46 . 

The results of our analysis is shown in Fig.l. The left-hand panel shows the contour 
plot in the mwiMP-^i plane for 'D=l of models Oi, i = 6,46,9,10 and 'D=1.7 for O 4 ; the 
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iiIdm (GeV) 


Figure 1. (a) Contour plot in the mwi m p-c^/<^ plane for the compatibility factor V defined in Eqs. 
(5. 1,5. 2). The constant value 'D=\ is shown for models Oi, i = 6,46, 9,10, while 22=1.7 is plotted for 
Oi- (b) For the same models Oi the minimum of V is plotted as a function of mwiMP when c”/cf 
is marginalized. 


right-hand panel shows the minimum of P as a function of mwiMP when rj is marginalized. 
A few features arise from both plots; i) when jcP —)-0 some intervals of mwiMP exist where 
the DAMA effect is compatible to constraints from other direct-detection experiments for 
generalized spin-dependent interactions {Oi, i = 6,46,9,10), while the compatibility factor 
for the standard spin-dependent interaction O 4 is always above P=l; ii) from Table 2 one 
can see that the pairs of models 0 ^-0^^ and Og-Ojo have the same momentum dependence, 
while they differ for the nuclear response function: since however — 2 Ws'/(g^) (at 

least for small q^) the compatibility factor P has very similar behaviours for the models 
in each pair ; ii) the tension between DAMA and other experiments is better alleviated in 
those models where momentum dependence is largest: in fact, models Oq and O 46 , which 
reach the best compatibility, depend on momentum through a factor compared to Og and 
Oio where the dependence is through q^ (see again Table 2); iv) for mwiMP > 30 GeV the 
compatibility factor rises steeply for models Oq, O 46 , Og and Oig, while no such feature is 
observed for the case O 4 . 

The last two properties can be understood in the following way. If mwiMP ^ 30 GeV, 
the WIMP signal in bubble chambers and droplet detectors is dominated by scatterings 
off fluorine. In particular, in this mwiMP range in order to deposit recoil energies above 
threshold scatterings off iodine in COUPP require Vmin values beyond the corresponding 
range for the DAMA signal, and we make the conservative assumption that in this case the 
halo function fj vanishes. Moreover scatterings off chlorine in SIMPLE are subdominant due 
to the suppressed nuclear response function. In this case the transferred momenta which 
explain the DAMA modulation effect for WIMP scatterings off sodium in DAMA are larger 
than the corresponding ones off fluorine in GOUPP, PIGASSO and PICO-2L. For instance, for 
?7iM5/MP=25 GeV one has 285 MeV^ < {q^AMA)'^ ^^70 MeV^, 275 MeV^ < (qcoupp)'^ ^^70 
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Figure 2. Measurements and bounds for the function fj defined in Eq. (3.4) for model O 4 (standard 
spin-dependent interaction) and c2=0. (a) mwiMP=W GeV; (b) mwiMP=^0 GeV. The / determi¬ 
nations for model O 7 would be rescaled by an approximately common factor ~ 3/2 with respect to 
the fj values shown in this figure. 

MeVE 60 MeV2 < (q^jcasso)^ ^470 MeVE 113 MeV^ < iq^jco- 2 L)^ ^V470 VleV^ . This 
implies that models where the expected detection rate depends on one additional factor 
(g^)”, n > 0 such as Oi, z = 6,46, 9,10 present a relative enhancement of the expected rate 
in DAMA compared to that in fluorine detectors, with a consequent relative loss of sensitivity 
for the latter. An exception to this argument is SIMPLE, where the interval for scatterings 
off fluorine has more overlap with that of DAMA, 283 MeV^ < qsiMPLE <470 MeV^, but 
which is overall less constraining than the other detectors due to the lower exposure (see 
Appendix B). On the other hand, for mwiMP > 30 GeV scatterings off iodine in COUPP 
become kinematically allowed, with values of the transferred momenta which are much larger 
than those related to fluorine and sodium. For instance, for mw/MP=35 GeV one has 1850 
MeV^ < {qaaiipp)"^ ^2350 MeV^. Clearly, this implies a strong enhancement of the expected 
signal in COUPP for interactions involving an additional dependence on 5 ^, with a consequent 
steep rise of the compatibility factor P, as observed in Fig. 1. This effect is not present for 
the standard spin-dependent interaction O 4 , whose expected rate has no q"^ dependence. 

In order to discuss in more detail the phenomenology in Figures 2, 3 and 4 the mea¬ 
surements and bounds for the function 77 obtained using Eq. (4.9) are plotted as a function 
of Vmin- In particular, in all these figures we take r* = 0. Fig. 2(a) shows the case of a 
standard spin-dependent interaction O 4 when mwiMP=^0 GeV, while Fig. 2(b) shows the 
same when mwiMP=^0 GeV. On the other hand. Figure 3(a) shows the case of model Oq for 
iPwiMP=W GeV while Fig. 3(b) shows the same when mwiMP=^^ GeV. Finally, in Figure 
2(a) and 2(b) model Og is shown for mwiMP=W GeV and mwiMP=^^ GeV, respectively. 

As discussed before, phenomenology of these models is practically degenerate with other 
scenarios. In particular, the ^ determinations for model Oj would be rescaled by an approx¬ 
imately common factor (~ 3/2) with respect to the fj shown in figure 2 for O 4 . In the same 
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way the phenomenology of O 46 is very similar to model Og, apart from an overall approxi¬ 
mate factor (~ 1/2) in the fj determinations, so it can be described by Figure 3. Moreover, 
the ^ determinations for model O 14 would be rescaled by an exact common factor of 2 with 
respect to the fj values shown in Figure 4 for Og; on the other hand, the fj determinations 
for model O 40 would be rescaled by an approximately common factor ~ 1/2. 

From all these plots we observe that germanium detectors and LUX imply relatively lax 
constraints, due to the suppression of their nuclear response functions when c/=0 (and, in 
the case of germanium, for the small natural abundance of the only isotope carrying 

spin). 

Moreover, in each plot the step-like cyan curve represents our conservative estimation 
fjl^^ivmin) of the minimal fji modulated halo function compatible to the DAM A signal. In 
particular, as explained before, we assume that fjl^^{vmin) vanishes outside the DAMA Vmin 
interval. The corresponding estimation fj’^^^iymin) > of the minimal fj function 

is obtained by requiring it to be a decreasing function of Vmin and is represented by the 
piecewise dotted red line. 

As explained before, we use the function fj^^^{vmin) to calculate the expected rate on 
droplet detectors and bubble chambers in the compatibility factor of Eq.(5.2). We decide to 
show the corresponding constraints by plotting for each experiment the function fj^^^{vmin) 
rescaled by a constant factor in such a way that the corresponding expected rate is equal to 
the corresponding most constraining upper bound on the number of events. As can be seen 
from all the figures, also in the most favorable situations {Oq ~ O 4 Q for mwiMP — 30 GeV) 
the bound is always relatively close (within one sigma) to the minimal DAMA signal. Since 
the DAMA estimations are for ?)i while the bounds are for 7) this implies that compatibility 
among them necessarily requires a modulation amplitude fraction in DAMA much larger 
then that predicted for a standard Isothermal Sphere (which is typically below 10%). In the 
same hgures the continuous cyan line and dotted red line represent and 

fjMaxweihan^y^^^-^ q£ g Maxwellian velocity distribution, where 

is normalized in such a way to minimize a y-square with the DAMA modulation amplitudes. 
Indeed, the red dotted line is well above existing constraints on fj in all cases. Moreover, notice 
that a Maxwellian for mwiMP=^0 GeV predicts a sizeable contribution to the event rate 
due to scatterings off iodine, so this explain the particularly poor fit with the fj estimations 
for DAMA, where dominance on sodium was assumed. 

7 Conclusions 

In the present paper we have used non-relativistic Effective Field Theory to classify the 
most general spin-dependent WIMP-nucleus interactions, and within this class of models we 
have discussed the viability of an interpretation of the DAMA modulation result in terms of a 
WIMP signal, using a halo-independent approach in which all dependencies from astrophysics 
are factorized in a single halo function. 

One of the main motivations of the spin-dependent scenario is the fact that the most 
stringent bounds on the interpretation of the DAMA effect in terms of WIMP-nuclei scat¬ 
terings arise today from detectors using xenon (XENONIOO, LUX) and germanium (GDMS) 
whose spin is mostly originated by an unpaired neutron, while both sodium and iodine in 
DAMA have an unpaired proton: if the WIMP effective coupling to neutrons is suppressed 
compared to that on protons this class of bounds can be evaded. In this case the most 
constraining remaining bounds arise from droplet detectors (SIMPLE, COUPP) and bubble 
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Figure 3. The same as in Figure 2 for model Oq. The fj determinations for model O 46 would be 
rescaled by an approximately common factor ~ 2 with respect to the fj values shown in this figure. 
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Figure 4. The same as in Figure 2 for model Og. The ^ determinations for model O 14 would be 
rescaled by a common factor 2 with respect to the 77 values shown in this figure. The fj determinations 
for model Oio would be rescaled by an approximately common factor ~ 1/2. 


Og DAMA = 

c|i=0 LUX / 

m^,„p=33 GeV SuperCDMS = 

XENONIO J 

CDMS-Light E 

CDMS-Ge 
SIMPLE E 

COUPP 

- PICASSO - I 

PICO-21 





10 -“ 

10-12 

10-13 
10-14 
10 “^ . 


10-16 
10 -“ 
1 n-iB 


a 


DAMA : 

LUX “ 

= 10 GeV SuperCDMS 

XENONIO m 

CDMS-Light^: 
CDMS-Ge 

_ 


SIMPLE 

COUPP 

PICASSO 

PICO-21 


cd- 


b E 


± '''■■4 (a) 


chambers (PICASSO, PICO-2L) , which all use nuclear targets (fluorine, chlorine and iodine) 
with an unpaired proton. 

From the phenomenological point of view we found that, although several relativistic 
Effective Field Theories can lead to a spin-dependent cross section, in some cases involving 
an explicit dependence of the scattering cross section on the WIMP incoming velocity (see 
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Table 2), three main scenarios can be singled out in the non-relativistic limit which (ap¬ 
proximately) encompass them all, and that only differ by their explicit dependence on the 
transferred momentum, n=0,l,2: using the notation of Eq.(2.2), they are O 4 , Og and 

Oe, respectively. 

In our quantitative analysis we pointed out that the requirement for a halo-independent 
analysis that one target dominates the expected rate is only needed to get estimations of the 
halo function from those experiments that observe an excess, but is not indispensable to 
get a conservative bound for those reporting a null result. We achieved this by adopting 
the procedure to find the minimal halo function compatible to the DAMA signal 

and then use it to calculate expected rates in other experiments, including droplet detectors 
and bubble chambers which contain several target nuclei. Using this approach we also con¬ 
cluded that scatterings on iodine can be assumed to be below threshold in DAMA as long as 
iT^wiMP ^ 60 GeV, allowing the factorization of the halo function for sodium and evading 
the KIMS bound (also using iodine) in the same WIMP mass range. 

In particular, we found that, for mwiMP ^ 30 GeV and with our assumptions on 
the WIMP signal in bubble chambers and droplet detectors is dominated by 
scatterings off fluorine. In this case models where the expected detection rate depends on 
one additional factor (g^)”, n > 0 show a relative enhancement of the expected rate in DAMA 
compared to that in fluorine detectors, with a consequent relative loss of sensitivity for the 
latter, because the transferred momenta which explain the DAMA modulation effect for 
WIMP scatterings off sodium in DAMA are larger than the corresponding ones off fluorine. 
In this way compatibility between DAMA and other constraints can be achieved for Og and, 
to a lesser extent, for Og, but not for the standard spin-dependent scenario O 4 . These 
conclusions are only valid for a WIMP velocity distribution in the halo of our Galaxy which 
departs from a Maxwellian. On the other hand, for mwiMP ^ 30 GeV a strong tension 
between DAMA and COUPP arises for both Og and Og because scatterings off iodine in 
COUPP become kinematically allowed, with values of the transferred momenta which are 
much larger than those related to fluorine and sodium. 
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A WIMP response functions 


We collect here the WIMP particle-physics response functions introduced in Eq.(2.4) and 
adapted from [24]: 
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B Experimental inputs for the analysis 

In this Appendix we summarize the experimental inputs that we have used to evaluate the 
response function defined in Eq. (4.5) for each of the experiments included in our analysis. 
Whenever applicable we will follow the convention to indicate with Er the true recoil energy, 
with Eee the electron-equivalent energy {E^e = Q{Er)Er with Q the quenching factor) and 
with E' the visible energy, as introduced in Section 4. In the case of bolometric measurements 
we assume Q = 1. With the exceptions of LUX we model the energy resolution with a 
Gaussian and we indicate the corresponding variance. 

DAM A We have taken the modulation amplitudes in 0.5 keVee bins from Fig .6 of Ref. [2] 
(already normalized to counts/day/kg/keV for a total exposure of 1.17 ton yr), adopting the 
signal region 2 keVee < E' < 4 keVee. We have adopted the value (5wa=0.3 for the quenching 
factor of sodium. 

LUX In the case of LUX we have assumed zero WIMP candidate events in the range 
2 PE< Si <30 PE in the lower half of the signal band, as shown in Fig. 4 of Ref. [3] for 
the primary scintillation signal Si (directly in Photo Electrons, PE) for an exposure of 85.3 
days and a hducial volume of 118 kg of xenon. Following Ref. [30] (see Eqs. (14-15)) we 
have modeled the detector’s response with a Poissonian fluctuation of the Si scintillation 
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photoelectrons combined with a Gaussian resolution crpMT=0.5 PE for the photomultiplier 
so that the response functions dehned in Eq.(4.5) are modified into: 


7^ 




8NT^^T^ Y [q‘^{Vmin)] 

T k rr' 


^Tk [yjVmin)] 
2jr + 1 



oo 

dSi Y Gauss (S'lln, y/nopMT) Poiss [n, v{Er)] icuts{Si), 

n=l 


m = 0,1. (B.l) 


In the equation above Poiss{n,X) = A”'/n! exp(—A), while ^cuts represents the combination 
of a 50% acceptance combined with the quality cut efficiency (taken from Eig. 9 of Ref. [3]). 
Moreover the expected number of PE for a given recoil energy Er is given by: 


k'iER) = Er X Leff{ER) X Ly, (B.2) 

with Ly=8.8 PE. We have taken L(,ff{Ef^ from [31] where it is calculated including the effect 
of the electric field, and assumed to vanish for S*! < 3 PE. 

SuperCDMS We include the low-energy analysis of SuperCDMS [10] with a germa¬ 
nium target in the energy range 1.6 keVnr< Er < 10 keVnr with a total exposition of 577 kg 
day and 11 observed WIMP candidates. The energy resolution is given by (TcDMS-Si{P') = 
y^O .2932 -I- 0.0562(E'/keVnr) in keVnr[32] 

XENONIO The analysis of XENONIO makes use of the secondary ionization signal 
S 2 only, with an exposition of 12.5 day and a fiducial mass of 1.2 kg. We take the scale of 
the recoil energy Er and the recorded event spectrum in the energy range 1.4 keVnr< Er < 
10 keVnr directly from Eig. 2 of Ref. [5]. The energy resolution is given by: (Txenonio = 
Er/ ^ERQy{ER) where Qy{ER) is the electron yield that we calculate with the same choice 
of parameters as in Eig. 1 of [5]. 

CDMSlite CDMSlite[9] analyzes the very low range 0.170 keVee< Eg^ <7 keVee for 
the electron-equivalent energy using a fiducial mass of 0.6 kg of germanium and an exposition 
of 10.3 days. We take the spectrum from Eig. 1 of Ref. [9]. We adopt the same quenching 
factor that we use for CoGeNT, an energy resolution acDMSUte =14 eV and the efficiency 
Ut =0.985 [9]. 

CDMS-Ge We consider the data from detector T1Z5 in the range 2 keVnr< Er < 100 
keVnr available in digital format from [8] with a raw exposure of 35 kg day on germanium 
target. The energy resolution is the same as in SuperCDMS, while the efficiency is taken 
from Eig.l of Ref. [8]. 

SIMPLE The SIMPLE experiment [11] uses superheated liquid droplets homogeneously 
distributed in a gel to search for transitions to the gas phase produced by WIMP scatterings. 
The nuclear targets are made of C 2 CIE^ (for the nuclear response function of chlorine, which 
is not available from [24, 27], we have used a simple estimation outlined in Appendix C). 
SIMPLE is a threshold detector, only sensitive to the minimal deposited energy Eth required 
to trigger the nucleations, and with Eth controlled by the pressure of the liquid. The prob¬ 
ability that an energy deposition Er on the target nucleus T nucleates a droplet is given 
by: 


Er - Eth 

Eth 


Pt{Er) = 1 - exp 


— CtT 


(B.3) 









Eth (keV) 

exposure (kg day) 

expected background (events) 

measured events 

95% upper bound 

7.8 

55.8 

0.8 

2 

5.92 

11 

70 

0.7 

3 

8.26 

15.5 

311.7 

3 

8 

12.29 


Table 3. For each operating threshold used in COUPP we provide the corresponding exposure, 
expected background, number of measured events and 95% C.L. upper bound obtained with the 
Feldman-Cousin method [33] used in our analysis 


Eth (keV) 

95% upper bound 

1.7 

1.1 

2.9 

1.5 

4.1 

11 

5.8 

9 

6.9 

1.3 

16.3 

3.1 

39 

1.5 

55 

6 


Table 4. 95% C.L. upper bounds (extracted from Fig. 5 of Ref.[13]) for each operating threshold 
used in PICASSO. 


where aT is determined by fitting calibrations with neutron sources. With an exposure of 
6.71 kg day and Eth=7.8 SIMPLE observed 1 event, consistent to an expected background 
of 2.2. This can be converted to an upper bound of 3.16 events using the Feldman-Cousin 
method [33]. We use aF=Oic=3.G. 

COUPP The COUPP experiment[12] searches for WIMPs using nucleations in a bub¬ 
ble chamber and is also a threshold detector. In the case of COUPP the target material 
is CF^I. For each operating threshold used in COUPP the corresponding exposure, ex¬ 
pected background, number of measured events and 95% C.L. upper bound obtained with 
the Feldman-Cousin method [33] used in our analysis are summarized in Table 3. We adopt 
the nucleation probability (B.3) with ai?=Q!c=0.15, while for iodine we assume P/=l, cor¬ 
responding to a/ —)• oo in (B.3). 

PICASSO The Picasso experiment [13] is a bubble chamber using CsFg, operated with 
eight energy thresholds. For each of the latter we provide the corresponding upper bound 
on the number of events (normalized to events/kg/day) in Table 4 (extracted from Fig. 5 of 
Ref. [13]). We use the nucleation probability of Eq.(B.3) with ac=aF=^- 

PICO-2L The PICO-2L collaboration operated a CsTg bubble chamber experiment 
with four energy thresholds. For each of them we provide the corresponding exposure, number 
of measured events and 95% C.L. upper bound (conservatively assuming zero background) 
used in our analysis in Table 5. In particular we conservatively chose to use the raw data 
without the subtraction adopted in [14] which makes use of time correlations among measured 
events. We adopt the nucleation probability (B.3) with aF=ac=0-15. 
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Eth (keV) 

exposure (kg day) 

measured events 

95% upper bound 

3.2 

74.8 

9 

16.77 

4.4 

16.8 

0 

3.09 

6.1 

82.2 

3 

8.25 

8.1 

37.8 

0 

3.09 


Table 5. For each operating threshold used in PICO-2L we provide the corresponding exposure, 
number of measured events and 95% C.L. upper bound (assuming zero background) used in our 
analysis. 


C Nuclear response functions for chlorine 

In the case of chlorine a shell model calculation for the nuclear response functions W^,' and 
is not available, so we assume = 2W^Ti and use a Gaussian approximation for 

the dependence. In particular, combining the usual spin-dependent scaling law written as 
[28]: 


cfn\ — ^ + l)(iT + 1) ^ c ^ I „ ^ o ^ ^2 

o (0) — . (dp ^ Op ^ On ^ j 

TT JT 


with the Gaussian form factor [34]: 


5(r 


5(0) 


(where At is the target nucleus mass number) the identities of Eq.(2.7) imply: 


(C.l) 


1 = ^ ^ + 2.68 - 0.78^ {A^^ - 3.8)2 + o.2^ fm, (C.2) 


('? ) — 
w^r'iq^) = 


8 (2jT + l)(jT + 1) 


37r 


jT 


< >< S^' > 


4 (2jT + l)(jr + 1) 


37r 


JT 


< 5^ >< 5^ > e 


IA 


(C.3) 


with < >= (< 5p > + < 5n >)/2 and < 5^ >= {< Sp > — < Sn >)/2. For both ^®G1 

and ^'^Cl we adopt < Sp >=-0.051 and < Sp >=-0.0088 [35]. 
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